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ABSTRACT 

We show that a relativistic gamma-ray burst (GRB) jet can potentially pierce the envelope of very 
massive first generation star (Population III; Pop III) by using the stellar density profile to estimate 
both the jet luminosity (via accretion) and its penetrability. The jet breakout is possible even if the 
Pop III star has a supergiant hydrogen envelope without mass loss, thanks to the long-lived powerful 
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accretion of the envelope itself. While the Pop III GRB is estimated to be energetic E 1 



10 



erg, the supergiant envelope hides the initial bright phase into the cocoon component, leading to a 
GRB with a long duration ~ 1000(1 + z) sec and an ordinary isotropic luminosity ~ 10 52 erg s _1 
(~ 10 -9 erg cm" 2 s" 1 at redshift z ~ 20). The neutrino-annihilation is not effective for Pop III 
GRBs because of a low central temperature, while the magnetic mechanism is viable. We also derive 
analytic estimates of the breakout conditions, which are applicable to various progenitor models. The 
GRB luminosity and duration are found to be very sensitive to the core and envelope mass, providing 
possible probes of the first luminous objects at the end of the high redshift dark ages. 

Subject headings: Cosmology: dark ages, reionization, first stars - Gamma-ray burst: general — Stars: 
general 



1. INTRODUCTION 

The ancient era of the first generation stars (Pop- 
ulation III; Pop III) - the end of the dark age - 
is still an unexplored frontier in the modern cos- 
mology dBarkana fe Loebll200H iBromm fc larsonl 12004 
ICiardi fc Ferraral 120051) . The first star formation from 
the metal free gas has a crucial influence on the sub- 
sequent cosmic evolution by producing ionizing photons 
and heavy elements. Although the theoretical studies 
are recently developed by the numerical simulations, the 
faint Pop III objects are difficult to observe even with 
the future technology. 

Gamma-Ray Bursts (GRBs) are potentially powerful 
probes of the Pop III era. In fact, GR B 090423 got 
the h i ghest redshift z = 8.2 ever seen (e.g.. iTanvir et al.l 
120091 : ISalvaterra et al.1 [20091 : IChandra et all l2010T >. be^ 
yond any qu asars or galaxies and previous GRB 080913 
at z = 6 .7 dGreiner et all 120091) and GR B 050904 at 
z = 6.3 (jKawai et all 120061 : iTotani et all 120061 ). The 
GRBs are presumed to manifest the gravitational col- 
lapse of a massive star - a collapsar - to a black hole 
with an accretion disk, launchi ng a collimated outflow 
(jet) with a relativistic speed (jMacFadven fe Wooslevl 
fl999h . The massive stars quickly die within the Pop III 
era. The GRBs, the most luminous objects in the Uni- 
yerse, are detectable in p rinciple out to redshifts z ~ 100 
(jlamb fc Reichartl 120001) , while their afterglows are ob- 



servab le up to z ~ 30 (ICiardi fc~loebl 120001: iGou et all 
120041: lloka fc Meszarosl 120051: iToma et all j2010D. As 
demonstrated in GRB 50904 by Subaru £k awai et all 
120061 : ITotani et all 120061 ) . the GRBs can probe the in- 
terstellar neutral fraction with the lyman a red damp- 
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ing wing (|Miralda-Escudelll998D. the metal enrichment, 
and the star formation rate iTotairj 119971 iKistler et all 
120091) . In the futur e, we may also inv estigate the reioniza- 
tion history (llokal 120031 : IInouell2004D . the molecular his- 
tory (llnoue et al.ll2007|), the equation of sta te of the Uni- 
verse (|Schaeferl 120071: lYonetoku et all 120041) and the ex- 
tragalactic backg round light (|Ohll200 li llnoue et al.ll2010b 
lAbdo et alJl2?nnt) . 

The first stars are p redicted to be predominantly 
very massive > 100M Q (|Abel et all 120021: IBromm et all 
2002). The mass scale is roughly given by the Jeans 
mass (or Bonnor- Ebert mass) when the isothermal- 
ity breaks (i.e., only through the cooling function 
and hence seems robust (but see also iTurk et~aTll200 
iClark et all l2010t ). The central part collapses first 
to a tiny (~ 0.01M Q ) protostar, followed by the 
rapid accretion of the surrounding matter to form a 
massive first star (jOmukai fc Pallal Il2003t lYoshida et all 
I2008D . The stars with 140-260M Q are expected to 
undergo the pair-instability supernovae without leav- 
ing any compact remnant behind, while those above 
~ 26OM would collapse to a massive (~ 100M Q ) 
black hole with an accreti on disk, potential l y leading to 
scaled-up collapsar GRB s (jFrver et al 2001; Heeer et all 



2003; 



2010; 



Suwa et all l2007a]lbl 120091 : iK omissarov fc Barkov 
Meszaros fc Reesl 120101) The Pop III GRB rate 
would be rare ~ 0.1-10 yr -1 but within rea ch (e.g., 
IBromm fc ioebl 120061 : iNaoz fc Brombergl l2007t ) . These 
GRBs also marks the formation of the first black holes, 
which may grow to the supermassive black h oles (BHs) 
via merger or accretion (|Madau fc Reesll200lD . 

However, the zero-metal stars could have little mass 
loss by the line driven wind (jKudritzkil 120021) . and 
thereby have a large (R* ~ 10 13 cm) hydrogen envelope 
at the end of life (red supergiant (RSG) phase). Es- 
pecially for Pop III stars, the mass accretion continues 
during the main sequence phase, so that the chemically 
homogeneous evolutio n induced by rapid rotation (e.g., 
lYoon fc langen 120051 ) might not work (jOhkubo et all 
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l2009f h Their extended envelopes may suppress the emer- 
gence of relativistic jets out of the ir surface even if such 
jets were produced ( Matznerll2003l ). The observed burst 
duration T ~ 100 s. providing an estimate for the life- 
time of the central engine, suggests that the jet can only 
travel a distance of ~ cT ~ 10 12 cm before being slowed 
down to a nonrelativistic speed. This picture is also sup- 
ported by the nondetections of GRBs associated with 
type II supernovae. Nevertheless, this may not apply 
to the Pop III GRBs because the massive stellar accre- 
tion could enhance the jet luminosity and duration and 
therefore enable the jet to break out the first stars. 

In this paper, we discuss the jet propagation in the first 
stars using t he state-of-the- a rt Po p III stellar structure 
calculated bv lOhkubo et al.l ([2009D (§[2j to estimate the 
jet luminosity via accretion (§ 12131) and to predict the ob- 
servational main characters of the Pop III GRBs, such as 
energy and duration. We adopt the analytical approach 
that reproduces the previous numerical simulations to see 
the dependences on the yet uncertain stellar structure 
(§ El f° r analytical estimates) and to avoid simulations 
over many digits. We determine the jet head speed that 
is decelerated by the shock with the stellar matter (§|4|. 
The shocked matter is wasted as a cocoon surrounding 
the jet before the jet bre akout (§[H), like in the co ntext 
of active galactic nuclei (jBegelman fc Ci offi 1989). We 
treat both the jet luminosity and its penetrability with 
the same stellar structure consistently for the first time. 

2. PROGENITOR STRUCTURE 

In this paper we employ three representative progen- 
itors; Pop III star, Walf-Rayet (WR) star, and RSG. 
These stars correspond to progenitors of Pop III GRBs, 
ordinary GRBs, and core-collapse supernovae without 
GRBs, respectively. The WR stars have no hydro- 
gen envelope, whi ch is a preferre d condition for a suc- 
cessful jet break (|Matznerl 120031 ) and consistent with 
the observational eviden ce of GRB-SN Ibc association 
(IWooslev fc Bloomll20M) . 

The density profiles of investigated models are shown 
in Fig. [T] Red line shows the de nsity profile of Pop II I 
star with 915 M , model Y-l of lOhkubo et al.1 (|2009h . 
Blue in dicates the GRB progeni tor with 16 M®, model 
16TI of IWooslev fc Hegerl ([2006). Green line represents 
the progenitor o f ordinary core-collap se supernovae with 
15 M , 815.0 of lWooslev et al.l ([20021 ). The density pro- 
files are roughly divided into two parts: core and en- 
velope. The GRB progenitor (WR star) does not have 
hydrogen envelope, while Pop III and RSG keep their en- 
velope so that these stars experience the envelope expan- 
sion triggered by core shrinkage after the main sequence. 

Because the exact stellar surface is difficult to calculate 
for the simulation of stellar evolution (K. Nomoto, pri- 
vate communication), we numerically solve the equation 
of hydrostatic equilibrium, 
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for the outermost layer of stars, where P is the pressure, 
r is the radius from the center of the star, G is the grav- 
itational constant, M r is the mass inside r, and p is the 
density, respectively. We employ the polytropic equation 
of state, P = Kp 1 , where K is the coefficient depending 
on the microphysics and 7 is the adiabatic index. Here, 
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Fig. 1. — Density profiles of investigated models. Red, blue, 
and green lines correspond to Pop III star (M = 915Mq), Walf- 
Rayet star (WR; GRB progenitor, M = 16Mq), and red supergiant 
(RSG; SN progenitor without GRB, M = 15M Q ), respectively. 
Pop III and RSG have a hydrogen envelope, which expands to a 
large radius, while WR has only a core. 

we use the constant value of K, fitting just outside the 
core. The surfaces of stars are determined by the point 
with P = 0. 

We can calculate the accretion rate, M, using these 
density profiles. The accretion timescale of matter at a 
radius r to fall to the center of the star is roughly equal 
to the free-fall timescale, 



tff 



GM r 



(2) 



Then we can evaluate the accretion rate at the cen- 
ter as M = dM r /dtff. Note that our estimation ne- 
glects the effect of rotati on (e.g., iKumar et al.l 120081 : 
iPerna feHYl acFadvcnl l2010h . However the rotation law 
inside the star is very uncertain. Even thought a ro- 
tationally supported disk is formed, the accretion time 
is roughly ~ a -1 = 10(a/0.1) _1 times tg, where 
a is the standard dimensionless viscosity parameter 
(|Kumar et al.l [20081) . In addition, the jet production 
mechanism is also unknown so that we introduce an ef- 
ficiency parameter to connect the (free-fall) mass accre- 
tion rate and jet luminosity, which will be normalized by 
the observed GRBs in the next section. This parameter 
would contain the information of both the rotation rate 
and the jet production efficiency. 

In Fig. [2l the mass accretion rates of investigated mod- 
els are shown. The origin of time in this figure is set at 
the time of BH mass (central mass) being 3M Q (3.4 sec 
after the onset of collapse for WR, for instance). The 
accretion rate should be related to the activity of the 
central engine, and that of Pop III stars is much larger 
than the other progenitors. Therefore, the GRBs of Pop 
III stars are expected to be more energetic than ordinary 
GRBs if Pop III stars could produce GRBs. However, it 
is nontrivial that the GRB jets can breakout the Pop 
III stars. In 21 we discuss the jet propagation and ca- 
pability of successful jet break. The colored regions in 
this figure show the hidden regions by the stellar interior 
where the jet propagates inside the star so that the high 
energy photons can not be observed. 

3. JET MODELS 
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Fig. 2. — Accretion rates as a function of time. Red, blue, and 
green lines show Pop III, WR, and RSG, respectively. Dotted re- 
gions represent that the jet propagates inside the star, while the 
solid regions correspond to the time after the jet breakout for the 
magnetic jet model. Solid lines give information of observables 
(e.g., duration and energetics of GRB). On the other hand, dot- 
ted regions show the hidden energy inside the star that goes into 
the nonrclativistic cocoon component. Th e gray dot-dashed line 
represents the analytic model in Eq. ( 1 1 3 1) ■ The black line shows 
j(3— 2n)/3 ^ £-0- 7 3 as a reference, where n ~ 2.6 is a parameter for 
the density profile (an effective polytropic index of the envelope; 
see 



In this study we employ the collapsar model, which 
is a widely accepted scenario for the central engine of 
long GRBs. In this scenario a black hole accompanied 
by the stellar collapse produces a relativistic jet, which 
is strongly suggested by observations. The greatest un- 
certainty in this scenario is the mechanism for convert- 
ing the accretion energy or BH rotation energy into the 
directed relativistic outflows. There are mainly two can- 
didates of the jet production in the vicinity of the central 
engine: neutrino annihilation and magnetohydrodynami- 
cal ( MHD) mechanisms includin g Blandford-Znajck pro- 
cess (jBlandford fe Znaieklll977t) . which converts the BH 
rotation energy into the Poynting flux jet via magnetic 
fields. Although t here are plentiful studies about these 
mechanisms (e.g., iPopham et al.|[l999t iDi Matteo et al.l 
12001 iProga et al.l 120031: iMcKinnevI 120061) . we have no 
concrete consensus for the available energy injection rate 
from the central engine into the jet. Therefore, we em- 
ploy two simple models for jet producing mechanisms. 
We assume that the jet injection luminosity can be writ- 
ten with functions of the accretion rate, M. The models 
used in this study is basically written in M or M 2 . The 
accretion-to-jet efficiency are given by the GRB observ- 
ables with the WR model in $5] More detailed expres- 
sions are following: 

(A) L ~ M model (MHD mechanism): A jet is driven 
by magnetic fields that are generated by accreting mat- 
ter 4 . In this case the jet injection luminosity is given 
by Lj = i)Mc 2 , where 77 is t he efficiency parameter. In 
iKomissarov fc Barkovl (|2010D . rj = 0.05/a/3, where f3 is 
the so-called plasma beta (/3 = SttP/B 2 with B being 
magnetic field), and we do not know the reliable values 
for both a and /3 in the collapsar scenario. Therefore we 

4 Although the existence of strong magnetic field in the Pop III 
stars is unclear, there are severa l studies on the magnetic field am- 
plification (e.g.,[Sur et al. 2 0101) . Here we assume that the strong 
magnetic field can be generated at the vicinity of the BH. 



parameterize these parameters with rj simultaneously. 

(B) L r~j M 2 model (neutrino-annihilation mecha- 
nism): A jet is driven by annihilation of neutrinos (yv — > 
e~e + ), which are copiously radiated by "hyperaccretion 
flow" (iMacFadven fc Woosievl 119991 ) . As for neutrino- 
annihilation process, the jet injection luminosity is writ- 
ten as Lj = CM 9/i M~^ /2 (jZalamea fc Beloborodovl 
I2010I) . where £ is the efficiency parameter including the 
information of accretion disk, e.g., the spectrum of emit- 
ted neutrinos and geometry of disk. 

4. PENETRATION OF STELLAR ENVELOPE 

In this section, we consider the propagation of the jet 
head in the progenitor star. If a relativistic jet (1^ ^> 1) 
strikes the stellar matter, two shocks are formed: a for- 
ward shock (FS) that accelerates the external material 
to a Lorcntz factor Th, and a reverse shock (RS) that de- 
celerates the head of the jet to Th- Balancing the energy 
density behind the FS (Pf) with that above the RS (P r ), 
one can obtain the Lorentz factor of the jet head. As for 
the ultra-relativistic case (T h ^> 1), Pf = ^T 2 pc 2 and 

P r = ^(^-) 2 njm p c 2 , where nj = L iso / '^r 2 T 2 c is the 
jet proper proton density with Li SO being the isotropic 
luminosity of a jet and m p is the proton mass, while 
for nonrelativistic case (Th «J 1), Pf = ^Y^-p(3 2 c 2 and 



Pr = 



±T 2 nj m p c 



where fit is the velocity of FS in unit 



of the speed of light, c. These equations leads follow- 

ing relations: ultra-relativist ic one, Th ~ L^r^ 1 / 2 p~ 1//4 
(M eszaros fc Waxmanl |2001[ ). and nonrelativistic one, 

Ph ^ Ll^r^p- 1 / 2 (|Waxman fc Meszarosl f2003h . Here 
we combine these equations empirically as follows: 

1/2 

hTl PS 18 ' 



/ L iso \' / r X- 1 
v 10 52 erg s- 1 J V10 12 cm/ 



P 



10 7 g cm 3 



-1/2 



(3) 



This approximation leads the same relation with 
IWaxman fc Meszarosl (|2003l) for the nonrelativistic c ase 
(r/j w 1) and agree with iMeszaros fc Waxmanl (|200ll ) to 
within 40% for the ultra-relativistic case (fih ~ !)■ The 
crossing time of the FS is also given by 



(4) 



As for the relativistic FS, the crossing time is much 
shorter than the light crossing time due to its large 
Lorentz factor (see IMeszaros fc Reesll200"ll) . 

Combining Eqs. (J3j> and ((4]), we obtain the necessary 
isotropic jet luminosity for the FS to reach the radius r 
as 



3 x 10 b 



10 12 cm 



10 7 g cm 3 

-2 



erg s 



(5) 



If the jet luminosity decreases slower than t~ 2 , the jet 
luminosity can achieve this value at the late phase. We 
can follow the evolution of the FS by equating Lj and 
Eq. ([5]) including the correction of the jet opening angle, 



(i.e., L, = L iso 6 2 /2). 
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We note that the accreting gas from the sur- 
rounding to the progenitor star is negligible for 
the jet breakout since the density of the accret- 
ing gas is low, p = M/Anr 2 v ~ 5 x 10~ 12 g 
cm" 3 (ill/lCT 2 M Q yr" 1 ) (r/10 13 cm)" 2 (u/10 8 cm s" 1 ) 

5. GRB AND COCOON 

In this section, we divide the energetics of the jet 
into two components: GRB emitter (relativistic com- 
ponent) and cocoon (nonrclativistic component). When 
the Lorentz factor of the FS, T/j, is smaller than 9J 1 , 
the sho cked material m ay escape sideways and form the 
cocoon (|Matznerll2003l) , which avoids the baryon loading 
problem. With this scenario, the injected energy before 
the shock breakout goes to the energy of cocoon and that 
after breakout goes to the GRB emitter. Therefore, we 
can calculate the energy budget of the GRB emitter and 
cocoon after the determination of the jet breakout time, 
t\y. We define tb as the maximum time obtained by Eq. 

First of all, we determine the accretion-to-jet conver- 
sion efficiency (depending on the mechanism) using the 
ordinary GRB progenitor (WR) to make the total energy 
of GRB emitter E to t = 10 52 erg 5 . In this estimation we 
assume that the half opening angle of the jet 9j = 5°. 
A successful GRB requires following two conditions; i) 
The jet head reaches the stellar surface, ii) The velocity 
of the jet head, shou l d be larger than that of the 
cocoon, (3 C (|Matznerll2003t iToma et al.ll2007t) . 

As for L ~ M model, the results of the WR case are 
Lj = 1.1 X 10 51 (Af/M Q s" 1 ) erg s" 1 , i.e., 

rj = jf- 2 « 6.2 X 10- 4 , (6) 

and t b = 4.7 s. For L ~ M 2 model, Lj = 76 x 10 51 
(M/Mq s- 1 ) 9 / 4 (M BH /M q )- 3 / 2 erg s" 1 and t b = 2.8 s 6 . 
We estimate the expected duration of the burst with the 
period during which 90 percent of the burst's energy is 
emitted, Tqq. Both models reproduce the typical dura- 
tion of burst of ~ 10 s (see Table [J). The energy of 
cocoon (injected energy before the shock breakout) is 
smaller than that of the GRB emitter. The isotropic 
kinetic energy of the GRB emitter is ~ 10 54 erg. 

Next, we apply the above scheme and jet luminosity 
(e.g., the same r\ and £) to the RSG (progenitor of su- 
pernovae without GRBs) and find that RSG can not 
produce GRB. This is because the jet head is slower 
than the cocoon. As for WR with L ~ M model, 
(i h ~ R*/(ct b ) ~ 0.3 and (3 C ~ ^/E c /(Mc 2 ) ~ 0.01, 
where E c ~ 2 x 10 51 erg is the energy of the cocoon 
(see Table 1) and M ~ IOM0 is the stellar mass, hence 
Ph > Pc On the other hand, /3 h ~ 0.007 and C ~ 0.01, 
i- e - fih < /3c, for the RSG. Thus the morphology of the 
shock wave is almost spherical and the jet cannot break 

5 Note that -Etot is not the total energy of gamma rays because 
there must be a conversion from the jet kinetic energy to gamma 
rays. Though the efficiency of conversion is unclear, it is typically 
the order of 0.1. Therefore, we employ -Etot = 10 52 erg that could 
lead the true gamma-ray energy of GRB ~ 10 51 erg. 

6 This luminosity sh o ws a similar value with Eq. (22) of 
Zalamca & Bcloborodov (2010) because if Mbh = 3Mq, Lj ss 
16 X 10 51 erg s" 1 . 



out the stellar surface with a small opening angle. In 
addition, the FS cannot reach the stellar surface with 
L ~ AI 2 model for the RSG. Therefore our scheme is 
consistent with observations of GRB-SN Ibc connection. 

Finally, we calculate the evolution of the jet head for 
the case of Pop III star (see Table [1]). We find that 
the L ~ M 2 model docs not produce GRB because the 
FS stalls inside the envelope due to rapidly decreasing 
jet luminosity (so-called "failed GRB"). On the other 
hand, the L ~ M model can supply enough energy for a 
jet to penetrate the envelope and produce a GRB. Since 
@h ~ 0.4 and j3 c ~ 0.08 in this model, the relativistic jet 
can penetrate the stellar envelope with a small opening 
angle and produce a successful GRB. The total energy 
of the GRB jet (injected energy after breakout) is ~ 45 
times larger than the ordinary GRB and the duration is 
much longer (Tgo ~ 1000 s). In addition, we estimate 
the minimum rj for the successful breakout, which is 77 w 
3.4 x 10~ 5 . This is 20 times smaller than that of the 
normal GRB. In this case (3 h - 0.03 and j3 c - 0.008. 
Below this value, the jet head cannot reach the stellar 
surface. 

The accretion of the envelope (not core) is very im- 
portant for the Pop III GRB. Although the envelope is 
mildly bounded by the gravitational potential (because 
7 w 1.38 ~ 4/3) and easily escapes when it is heated by 
the shock, the timescale of cocoon passing in the enve- 
lope (~ J?*/(c/3 c ) ~ 3000 s) is longer than tf,. Therefore, 
the envelope accretion can last till the jet breakout and 
our conclusion about the penetrability of the relativistic 
jet is not changed. 

It should be noted that the opening angle of the jet 
could not be constant during the propagation phase. Due 
to the additional collimation by the gas pre ssure, 8 j be- 
come s smaller as the jet propag ates (e.g., IZhang et al.l 
120031 : iMizuta et all I2006L l2010t ). The smaller 6j leads 
the larger L iso (= 2Lj/9 2 ) so that our constant 6j is a 
conservative assumption for the jet breakout. 

6. ANALYTICAL DEPENDENCES ON PARAMETERS 

The hydrogen envelope could be reduced by the mass 
loss, which is one of the most uncertain processes in the 
stellar evolution. Even in the zero-metal stars, the syn- 
thesized heavy element could be dredged up to the sur- 
face, and might induce the line driven wind. The stellar 
luminosity could also exceed the Eddington luminosity of 
the envelope. The stellar pulsation might also blow away 
the envelope dynamically. So we analytically estimate 
the dependence on the envelope mass in the following. 

The density profile of the envelope can be written as 

where n is a constant (jMatzner fc McKeelll999l ). This 
profile is exact if the enclosed mass is constant (i.e., the 
envelope mass is negligible compared with the core mass) 
and the equation of state is polytropic, in which case 
n = (7 — 1) _1 is a polytropic index. We have n = 3/2 for 
efficiently convective envelopes since the adiabatic index 
is 7 = 5/3 for the ideal monoatomic gas, and n = 3 for 
radiative envelopes of constant opacity k since P oc p 4 / 3 
is derived from the relations, P 7 = (L/LEdd)P, P pT, 
P 7 cx T 4 , with a constant luminosity- mass ratio L/M 
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where LEdd = 4nrGMc/ n > L is the Eddington luminos- 
ity. We can fit well the Pop III envelope in Fig. Q] with 
n w 2.6. 

Using the profile in equation (JT)), we can estimate the 
envelope mass as 



p(r)4:Trr 2 dr oc 



Pi 



ir 



R" Ri- 



ft 



where p c = p(R c ) ~ p\{R* / R c ) n is the envelope density 
just above the core and R c is the core radius. The core 
density is higher than p c to proceed with the nuclear 
burning. Since the density enhancement is determined 
by the difference of the ignition temperature, which is 
not sensitive to other parameters, we assume that the 
core mass is given by M c oc p c Rc, so that 



M C R?~ Rl 



(9) 



Then the stellar radius can be written as a function of 
the core radius, the core mass, and the envelope mass as 
follows: 



( ^"env 



V 500 Ma 



-2.5 



2.5 



(10) 



where we use n = 2.6 (see Appendix for n dependences). 
Therefore, the stellar radius has a strong dependence on 
the envelope mass. If the envelope mass is smaller than 
~ 50Mq, the stellar radius is almost the core radius, 
R c ~ 10 10 cm. 

Next, we consider the jet breakout time. Since the 
accretion time of the core (r < 10 10 cm) is ~ 4 s, the en- 
velope accretion is important for the successful breakout 
if tb is longer than this timescale. Using iff ~ -J r 3 /GM r , 
we can evaluate the envelope accretion rate as 

M = MIr/dr oc Pl Mi 3 -")/ 3 <i (3 - 2n)/3 (11) 
dtfi jdr 

with the approximation of M r — M c = 
f R p{r')A-Kr' 2 dr' < M c 400M Q , which is valid for 

r < 10 12 cm (corresponding to ts < 3000 s). Combining 
with Eqs. ©, (dOjl, (dU and Lj = rjAIc 2 = L iso 9 2 /2, 
the breakout time is given by 



t b ~ 700 s 



V 



-0.79 / Q 



M, 



lO- 3 / 

2.!) 



1.6 



R c 



400M, 



V 500M, 







10 10 cm 

2.8 



1.1 



(12) 



where we use n = 2.6 (see Appendix for n dependences) 
and put p = p\ and r = i?„ in Eq. ([5]). If 77 is very small 
(< 10~ 5 ) and tb is longer than the free-fall timescale 
of the outermost part of star (~ 10 5 s), the jet cannot 
penetrate the star. When 77 sa 10 -5 , fth ~ R*/h ~ 
0.003 and f3 c ~ ^/rj ~ 0.003 so that the shock wave also 
propagates almost spherically. Similarly, the RSG case in 
Fig. [5] also has too long breakout time since M c ~ 4M Q 



and M e 



The jet luminosity after the breakout (t > tb) is given 



by 



,(t)~5 x 10 52 ( 



V 



Mr 



10- 3 
1.1 



Rr 



10 1 



/ t 



4ooM y V 700 s 



-0.73 



erg s 



(13) 



where we use n = 2.6 (see Appendix for n dependences) 
and M c oc p c R^ — Pi(R*/ 'R c ) n Rl- Interestingly the de- 
pendence of Li so (r;) on the envelope mass is only through 
the time, t. The relativistic jet emitted after the break- 
out will produce GRB prompt emission and the jet lumi- 
nosity decreases with time as /j~( 2 "~ 3 )/ 3 ~ £~ - 73 in this 
case. Fig. [5] shows this dependence with the gray dot- 
dahsed line, which well reproduces the numerical result 
(red line) for 10 s < t < 3000 s. The parameter de- 
pendences just after the breakout (most luminous time; 
t = tb) are given by Eqs. (fT2j) and (fT3|) as 

-3.2 / n \ -1.2 



L iso (t = t b )~5x 10 



52 



(lO- 3 ) 



1.6 Q 



R, 



( M c 



V400M G 



3.2 



/ M c 



V500A/ G 



10 10 cm 



erg s 



(14) 



T, 



■ in ) 



where we use n = 2.6 (see also Appendix). 

The active duration of the central engine (as % 
is determined by the accretion timescale at r ~ 10 12 cm, 
where the density gradient gets larger because p(r) oc 
(R*/r — 1)". The region of r < 10 12 cm has p(r) oc 
{R*/r) n so that Eqs. ([12]) and are valid, while they 
are not appropriate for r > 10 12 cm (M r > M c +0AM env ) 
due to the existence of the stellar surface, which leads the 
fast decrease of the accretion rate (compare the red and 
gray lines in Fig. [5]). So, the duration can be calculated 
by t ff at r — 0.1R* with Eq. (fTDJ as 



< ff (r = 0.LR*) ~ 3000 



10 



1(1 



1.5 



M c 



AOOMr. 



-3.8 



3.8 



M c + 0.4M e 



-0.5 



(15) 



500MqJ V 6OOM 

where we use n = 2.6 (see also Appendix). This is con- 
sistent with the result tb + T90 ~ 2200 s in the previous 
section. At this time, the isotropic luminosity in Eq. ([Ti 
is given by 

L iso [t = t s (r = 0.LR*)] ~ 2 x 10 52 (%L 



Rr. 



10 



10 



-1.5 



M c + 0AM e 



V400M Q 
0.37 



3.0 



V500A/ Q 



-2.8 



600A/, 







erg s 



(16) 



where we use n = 2.6 (see also Appendix). We note 
that the observables in Eqs. (fbi]) . (fl~5j) and ([T5]) carry 
information of the stellar structure. 

7. SUMMARY AND DISCUSSION 

We have investigated the jet propagation in the very 
massive Population III stars assuming the accretion-to- 
jet conversion efficiency of the observed normal GRBs. 
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We find that the jet can potentially break out the stellar 
surface even if the Pop III star has a massive hydrogen 
envelope thanks to the long-lasting accretion of the enve- 
lope itself. Even if the accretion-to-jet conversion is less 
efficient than the ordinary GRBs by a factor of ~ 20, 
the jet head can penetrate the stellar envelope and pro- 
duce GRBs. Although the total energy injected by the 
jet is as large as ~ 10 54 erg, more than half is hidden 
in the stellar interior and the energy injected before the 
breakout goes into the cocoon component. The large en- 
velope accretion can activate the central engine so that 
the duration of Pop III GRB is very long if the hydrogen 
envelope exists. As a result, the luminosity of Pop III 
GRB is modest, being comparable to that of ordinary 
GRBs. 

Considering the Pop III GRB at redshift z, the dura- 
tion in Eq. (| 15[) is 



Ti 



GRB 



T 90 (l + z) « 30000 s 



1 



20 



(17) 



which is much longer than the canonical duration of 
GRBs, ~ 20 s. The total isotropic-equivalent energy of 
Pop III GRB is 



1.2 x 10 55 (^1 
VO.I 



erg, 



(18) 



where e 7 is the conversion efficiency from the jet ki- 
netic energy to gamma rays (see Table []}. It should 
be noted that this value is comparable to the largest 
fi. iso ever observ ed, 9 x 10 54 erg for GRB 080916C 
(|Abdo et al. 2 0091). This value i s smal le r than the estima- 
tion of iKomissarov fc Barkovi (|20Tnh : M cszaros fe Reesl 
(|2010l ) because we consider the hidden (cocoon) compo- 
nent. Since the large isotropic energy is stretched over 
the long duration, the expected flux just after the break- 
out is not so bright, 



F 



10 9 erg cm 2 s , 



(19) 



where ri, is the luminosity distance, which is smaller 
than the Swift Burst Array Telescope (BAT) sensitivity, 
~ 10 -8 erg cm -2 s _1 . However, there must be a large 
variety of the luminosity as ordinary GRBs so that more 
luminous but rare events might be observable by BAT. 
Although the cocoon component has a large energy, the 
velocity is so low that it is also difficult to observe. If 
the cocoon component interacts with the dense wind or 



ambient medium, it might be observable. This is an in- 
teresting future work. 

The above discussions strongly depend on the enve- 
lope mass because the stellar radius is highly sensitive 
to the envelope mass. We derive analytical dependences 
on the model parameters in § |6l which are applicable to 
a wide variety of progenitor models. According to the 
analytical estimates, the smaller envelope leads to the 
shorter duration and the larger observable luminosity. 
If Af onv < 50M Q , R* ~ R c ~ 10 10 cm (sec Eq.(TO), 
tb ~ 2 s (obtained using numerical model as in fJSJ), 
and T 90 ~ t s ~ 4(M/400M o )- 1 / 2 (i? c /10 10 cm) 3 / 2 s; 
Both tb and T90 might be extended by a factor of a -1 
due the the rotation. This is much shorter than the 
GRBs from Pop III stars with massive envelopes. The 
mass accreting after the breakout is ~ 200M Q so that 
Ej.iso ~ 3 x 10 54 (e 7 /0.1) erg can be emitted by gamma 
rays after breakout. The luminosity just after the break- 
out is L- v 



10 erg s \ i.e. 



10 7 erg cm 2 s 1 , 



which is much brighter than the case with the massive 
envelope and observable with the Swift BAT, while the 
duration is very short (~ 2 s at the source frame). 

The matter entrainmcnt from the envelope is also cru- 
cial for th e fireball dynamics and the GRB spectra (e.g., 
IIokall20T0h . Since the envelope of the PopIII star is differ- 
ent from that of present-day stars, the GRB appearance 
is also likely distinct from the observed ones. This is also 
an interesting future work. 

Since we do not have any conclusive central engine sce- 
nario, we employ two popular mechanisms in this paper, 
that is, the magnetic model (L ~ M) and neutrino- 
annihilation model (L ~ M 2 ). The difference between 
these models is the dependence of the mass accretion 
rate. We find that the neutrino-annihilation model can- 
not penetrate the Pop III stellar e nvelope, so that n o 
GRB occurs, which is consist e nt withlFrver et alj d200ll) ; 
IKomissarov fe Barkovi (|2010f k iMeszaros fe Reed mm . 
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APPENDIX 

DEPENDENCE ON THE DENSITY INDEX N 
Here we explicitly show the dependences on n in Eqs. (fTTJ)) . (fl"2"|) . (fTS")) . (fl4"|) . (|15l) and ([TBI as follows, 

R* cx R C M~^ [(3 - n)Af cnv ] , 

6 3(4-n) „2_ 9 „ + 21 3(4-t0 

(3_„)(9_2„) 



M c 

2n—3 



[(3 - n)M e „v] (3-»>«>-*o 



£iso(*)<X770- 2 i?-( 3 -")M c 3 t~"^~, 



t s (r = 0.1R*)ocrIMc ] 



Rc 



M, 



[(3 - 7l)Af onv ] (3-n H 9-2„; 



L iso [t = ts(r = Q.lR^oz^e-'Rc 5 M c 



[(3 - n)M ew ]^(M c + 0.4M env )-i, 

~ [(3 - n)M cnv ] 5pSr (M c + 0.4M cnv ) ^ 



(Al) 
(A2) 
(A3) 

(A4) 
(A5) 
(A6) 
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TABLE 1 
Model Summary 



Model 


Final Mass 


Radius 


Mechanism 


*b 


Energy of GRB emitter 


Energy of Cocoon 


Tgo 


Eiso 


(Reference) 


[Mq] 


[10 11 cm] 




w 


[10 52 org] 


[10 52 erg] 


W 


[10 54 erg] 


WR 


14 


0.4 


MHD 


4.7 


1.0 


0.23 


49 


2.6 


(Wooslcy & Hcecr 2006) 






Neutrino 


2.8 


1.0 


0.42 


10 


2.6 


Pop III 


915 


90 


MHD 


690 


45 


57 


1500 


120 


(Ohkubo et al. 2009) 






Neutrino 












RSG 


13 


600 


MHD 


(5 x 10 5 ) 


(0.88) 


(0.16) 


8.4 x 10 6 ) 




(Wooslcy ct al. 2002) 






Neutrino 













Note. — The final masses of WR and RSG are not the same mass as the initial mass of progenitors (15Mq and I6M0) because of 
mass loss due to stellar wind, while Pop III stars do not have mass loss so that the final mass is 915Mq. t;, is the time of shock breakout 
from our calculations (see §(5j. The energy of GRB emitter is the kinetic energy of relativistic jet after the breakout. The energy of the 
cocoon is the injected energy before the breakout. E- lso is the isotropic energy of GRB emitter, corrected by the half opening angle of the 
jet, 8j = 5°. The parenthetic values for RSG are just for references because the cocoon velocity is larger than GRB emitter so that the 
GRB is not generated. 

We note that the total energy is proportional to 

tfr(r = QAR*)L iso [t = iff (r = OAR.)] oc rfiJ^M^ [(3 - n)M env ](M c + 0.4M cnv )-^. (A7) 
if n < 3 (i.e., L- iso (t) is shallower than tr 1 ). 
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